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^ ■ ABSTRACT 

^ \ The WFCAM Transit Survey (WTS) is a near-infrared transit survey running on the 

United Kingdom Infrared Telescope (UKIRT) , designed to discover planets around M 
dwarfs. The WTS acts as a poor-seeing backup programme for the telescope, and rep- 
resents the first dedicated wide-field near-infrared transit survey. Observations began 
K^ ■ in 2007 gathering J-band photometric observations in four (seasonal) fields. In this pa- 

ON ' per we present an analysis of the first of the WTS fields, covering an area of 1.6 square 

2^ ■ degrees. We describe the observing strategy of the WTS and the processing of the data 

to generate lightcurves. We describe the basic properties of our photometric data, and 
measure our sensitivity based on 950 observations. We show that the photometry 
^^ , reaches a precision of ^ 4 mmag for the brightest unsaturated stars in lightcurves 

^D ■ spanning almost 3 years. Optical (SDSS griz) and near-infrared (UKIRT ZYJHK) 

^^ \ photometry is used to classify the target sample of 4600 M dwarfs with J magnitudes 

in the range 11-17. Most have spectral-types in the range M0-M2. We conduct Monte 
Carlo transit injection and detection simulations for short period (<10 day) Jupiter- 
and Neptune-sized planets to characterize the sensitivity of the survey. We investigate 
the recovery rate as a function of period and magnitude for 4 hypothetical star-planet 
^ ; cases: MO-2-h Jupiter, M2-4-h Jupiter, MO-2-hNeptune, M2-4-hNeptune. We find that 

the WTS lightcurves are very sensitive to the presence of Jupiter-sized short-period 
transiting planets around M dwarfs. Hot Neptunes produce a much weaker signal and 
suffer a correspondingly smaller recovery fraction. Neptunes can only be reliably re- 
covered with the correct period around the rather small sample (~ 100) of the latest 
M dwarfs (M4-M9) in the WTS. The non-detection of a hot- Jupiter around an M 
dwarf by the WFCAM Transit Survey allows us to place an upper limit of 1.7-2.0 per 
cent (at 95 per cent confidence) on the planet occurrence rate. 

Key words: stars: planetary systems - stars: late-type ~ stars: statistics ~ infrared: 
stars - techniques: image processing 
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1 INTRODUCTION transiting exoplanet^ liave been discovered around F,G and 

K dwarfs, only ~ 70 sucli planets are known around the 
lower mass M stars. M dwarfs make interesting targets for 
many reasons, for example they provide better sensitivity to 

M dwarfs are the most numerous stars in our Galaxy 
l|Chabrieiil2003l ) . and until recently have remained relatively 
unexplored as exoplanet hosts. While several hundreds of ^ http://exoplanet.eu 



G. Kovdcs et al. 



smaller (rocky) planets in the habitable zone, and they also 
provide strong constraints on planetary formation theories. 

The most successful techniques so far used to discover 
exoplanets are the radial velocity (RV) and transit meth- 
ods. A planet around an early M2 dwarf can have ~4 times 
deeper transits and ~1.5 times higher RV amplitudes at the 
same orbital period than orbiting a Sun-like star. Both meth- 
ods also show increasing sensitivity towards smaller semi- 
major axis. Around solar-type stars, planets in short period 
orbits (days) are extremely hot (above lOOOK). If M dwarfs 
harbour planets in equally close-in orbits, they will proba- 
bly be more interesting due to their moderate equilibrium 
temperatures. In case of later (M4-9) type M dwarfs, these 
planets may fall in the star's habitable zone, and for rocky 
planets water can exist as a liquid, while Neptune-like plan- 
ets may have water vapour in their atmospheres. 

M dwarfs are also important for testing planet forma- 
tion theories. In the core accretion paradigm, dust parti- 
cles in the protoplanetary disk coagulate and form solid 
cores. The cores then begin to accumulate gas from the 
disk. When cores reach several Earth masses, gas accu- 
mu lation can speed up significantly . It w as first shown 
by iLaughlin. Bodenheimer fc Adamsl (|2004 ) that gas gi- 
ants cannot form easily around low mass stars this way. 
In the core accretion process gas giant formation is in- 
hibited by time scale differences. The protoplanetary disk 
around a low mass host dissipates before planetary cores 
are able to accrete their gas envelope by entering the run- 
away gas accretion ph ase. These planets remain 'failed 
cores'. Ilda fc LinI (|2005l ) found very similar results. Newer 
models refine this picture by introducing more detailed 
relations between stellar and disk properties. They con- 
sider migration driven by disk-planet interactions (type 
I-II) ([Ida fc LinI l 200g). evolution of snow-line location 
Kennedv fc Kenvo n 2008 ) and plan et-planet scattering 



around solar G ty pe dwarfs. RV studies determined a rate 
of 0.9-1.5 per cent JMarcv et al.ll2005l:[Cu"mming et al1l2008l : 



iMavor et al.ll201ll : IWright et al 



20121 '). while transit studies 



Thommes. Matsumura fc Rasid 120081 ). Their conclusions 



allow gas giant formation around low mass stars but the pre- 
dicted frequency of gi ants systematically decreas es towards 
lower primary masses. [Kennedy fc KenvonI (|2008l ) give rela- 
tive ratios for gas giant planeto ^ fraction of 1 per cent of 
low mass stars is predicted to have at least one giant planet, 
assuming that this ratio is 6 per cent for solar mass ones. 

Due to the failed core outcome, it is not surprising that 
Neptunes and rocky planets are predicted to b e common 
around low mass stars in core accretion models. Ilda fc LinI 
l|2005l ) predict the highest frequency around M dwarfs for a 
few Earth mass planets in close-in orbits (<0.05 AU). Later 
models also firmly support the existence of these smaller 
planets (|lda fc Linll2008l : iKennedv fc Kenvonll2008l '). 

Even if it is hard to form giant planets via core accre- 
tion around the lowest mass stars, gravitational instability 
models can produce gas giants on a very short timescale 
(~ 10^ yr). They predict that gas giants can form around 
low mass primaries as efficiently as around more massive 
ones assuming that the pro toplanetary disk is sufficiently 
mass ive to become unst able (|Bossll2006r ). 

(|Wright et al.ll2012l . and references therein) give a sum- 
mary of occurrence rates of hot giant planets (T < 10 days) 



^ Planet formation theories usually do not provide absolute num- 
bers because of free interaction coefficients in their formulae. 



have a systematicall y lower rate (roughly half of this value) 
at 0.3-0.5 per cent JGould et al.ll2006l : [Howard et al.ll2012l . 
H12 hereafter). 

For M dwar fs, recent RV studies support the paucity 
of giant planets (ICumming et al.ll2008l : Ijohnson et al.ll2007l . 



I2OI0I : iRodler et al.ll2012l ) though there are co nfirmed detec- 
tions b oth with short- and long-period orbits. Ijohnson et al.l 
I2OO7I) found 3 Jovian planets (with orbital periods of years) 
in a sample of 169 K and M dwarfs, in the California and 
Carnegie Planet Search data. The planetary occurrence rate 
for stars with M< O.7M0 in their survey is 1.8 per cent, 
which is significantly lower than the rate found around more 
massive hosts (4.2 per cent for Solar-mass stars, 8.9 per 
cent around higher mass subgiants). The positive correla- 
tion between planet occurrence and stellar mass remains 
after metallicity is taken into account, although at some- 
what lower significance. Gravitational microlensing pro- 
gra mmes also det e cted J upiter-like giants aro und M dwarfs 
(e.g. ICould et al.l (|2010l 'l: [Batista et all (120111 )) but statisti- 
cal studies seem to arrive at different conclusions. Microlens- 
ing surveys are more sensitive to longer period systems than 
RV (and particularly transit) surveys. iGould et al.l (|2010l ) 
analyzed 13 high magnification microlensing events with 6 
planet discoveries around low mass hosts (typically O.5Af0) 
and derived planet frequencies from a small but arguably 
unb iased sample. They com pared their planet frequencies to 
the ICumming et all l|2008h RV study. They found that af- 
ter rescaling with the snow-line distance to account for their 
lower stellar masses in the microlensing case, the planet fre- 
quency at high semi-major axes is consistent with the distri- 
bution from the RV study extrapolated to their long orbits. 
They found a planet fraction at semi-major axes beyond the 
snow line to be 8 times higher than at 0.3 AU. Considering 
that hot giants in RV discoveries (around solar type stars) 
are thought to have migrated large distances into their short 
orbits, they conclude that giant planets discovered at high 
semi-major axes around low mass hosts do not migrate very 
far. Their study suggests that rather than the formation 
characteristics, the migration of gas giants may be different 
in the low mass case. 

Focusing again on giant planets with short orbital pe- 
riods around M dwarfs, then as of writing, there is only 
one c onfirmed detection. The Kepler Mission JBorucki et al.l 
Il997l). includes a sample of 1086 low-mass targets (in Q2) 
( Borucki et al.l I2OIII . IH12| ) and one confirmed hot Jupiter 
(P=2.45 days) arou nd an early M dwarf host (KOI-254) 
(| Johnson et al.ll2012l ). Unfortunatel y, th is object was not in- 
cluded in the statistical analysis of lH12l study. 

For transit surveys, sufficiently bright M dwarfs make 
good targets because of their smaller stellar radius. A 
Neptune-like object in front of a smaller host can produce 
a similar photometric dip (~ 1 per cent) as a Jupiter- 
radius planet orbiting a Sun-like star. Ground based sur- 
veys are therefore potentially sensitive to planets around M 
dwarfs that could not be detected around earlier type stars 
with typical ground-based precision. On the other hand, M 
dwarfs provide a more demanding technical challenge. They 
are intrinsically faint and their spectral energy distribution 
peaks in the near infrared. Their faintness at optical wave- 
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lengths also makes spectroscopic follow-up observations diffi- 
cult. Measuring photometry in the infrared helps, but intro- 
duces a higher sky background level. It is also worth pointing 
out that the intrinsic variability of M dwarfs (flaring, spots) 
could further reduce the ease of discovering transiting sys- 
tems. 

One approach is to target a large number of M dwarfs, 
with a wide-field camera equipped with optical or near- 
infrared detectors. An alternative approach is to target in- 
dividual brighter M dwarfs in the optical deploying several 
small telescopes. Up till now, only 2 tran siting planets have 
been discovered around ( bright) M stars (|Gillon et al.ll2007l : 
ICharbonneau et al.ll2009l ') by ground based transit surveys. 

In this paper, we discuss the UKIRT (United King- 
dom Infra Red Telescope) WFCAM (Wide Field CAMera) 
Transit Survey (WTS), the first published wide-field near- 
infrared dedicated programme, searching for short period 
(<10 day) transiting systems around M dwarfs. The sur- 
vey was designed to monitor a large sample (~ 10, 000) of 
low-mass stars with precise photometry. In this paper we 
present an analysis of the first completed field in the WTS. 
We demonstrate that we can already put useful constraints 
on the hot Jupiter planet occurrence rate around M dwarf 
stars, and this is currently the strictest constraint available. 

The survey observing strategy is described in Section 
[21 In Section [3] a summary is given of the data processing 
pipeline, describing the generation of final clean lightcurves 
from raw exposures. In Section [4] we describe the M dwarf 
sample and discuss uncertainties in classification. We eval- 
uate the sensitivity of the survey in Section [5] using transit 
injection and detection Monte Carlo simulations. The lack 
of giant planets detected by the survey around M dwarfs to 
date is d iscuss ed in Section [6] Finally, in Section [T] we con- 
sider the IH12| sample and use it to place an upper limit on 
the frequency of hot Jupiters around M dwarfs based on the 
Kepler Q2 data rel ease. We discuss our WT S resu lts in the 
context of both the lBonfils et al.l (|201lh and lH12l studies. 



2 THE WFCAM TRANSIT SURVEY 

UKIRT is a 3.8m telescope, optimized for near-infrared ob- 
servations and operated in queue-scheduled mode. The Min- 
imum Schedulable Blocks are added to the queue to match 
the ambient conditions (seeing, sky-brightness, sky trans- 
parency etc). The WTS runs as a backup programme when 
observing conditions are not good enough (e.g. seeing > 1 
arcsec) for main surveys such as UKIDSS ([Lawrence et al.l 
I2OO7I ). 

The Wide Field Camera (WFCAM) comprises four 
Rockwell Hawaii-II PACE arrays, with 2kx2k pixels each 
covering 13.65 arcmin x 13.65 arcmin at a plate scale of 0.4 
arcsec/pixel. The detectors are placed in the four corners 
of a square with a separation of 12.83 arcmin between the 
chips. This pattern is called a pawprint. 

The WTS time series data are obtained in the J band 
(Aeff = 1220 nm), the fields were also observed once in all 
other WFCAM bands (Z,Y,H,K) at the begin ning of the sur- 
vey. T his photometric system is described in lHodgkin et al.l 
((20091). Each field, covering 1.6 square degrees, is made up of 
an 8 pawprint observation sequence with slightly overlapping 
regions at the edges of the pawprints (Fig.[T|). Observations 
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Figure 1. Observation pattern for the WTS. A field of 1.6 square 
degrees consists of 8 pawprints (a-h), each pawprint is built up 
from the simultaneous exposures of the four (numbered) detec- 
tors. The X-axis is Right Ascension, and increases to the left of 
the figure, while the Y-axis is Declination, and increases to the 
top. 
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Table 1. Summary of the WTS fields and their coverage as of 
27th May 2010. Stellar objects are morphologically identified by 
the photometric pipeline (see Sec l3.2l l. 



are carried out using 10s exposures in a jitter pattern of 9 
pointings. A complete field takes 16 min to complete form- 
ing the minimum cadence of the survey. Observing blocks 
typically comprise 2 or 4 repeats of the same field. 

Four distinct WTS survey fields were chosen to be rea- 
sonably close to the galactic plane to maximize stellar den- 
sity while keeping giant contamination and reddening at 
an acceptable level (see Section Q. The four fields are dis- 
tributed in right ascension at 03, 07, 17 and 19 hours to 
provide all year coverage (at least one field is usually visi- 
ble). At the time of the analysis presented in this paper, the 
19 hour field has approximately 950 epochs which is close to 
completion (the original proposal requested 1000 exposures 
for each field), the other three are less complete. A sumiuary 
of the key properties of the survey regions is shown in Table 

m 

The WTS has a lower priority than most of the main 
UKIRT programmes, thus observations are not distributed 
uniformly over time. For any given field, such as the 19 hour 
field (Fig. [2]), there are large gaps when the field is not visi- 
ble, as well as variations between and within seasons. 



3 PIPELINE OVERVIEW 

The WTS uses list-driven aperture photometry on processed 
images to construct lightcurves from the stacked data frames 
at the pawprint level. We describe this procedure briefly in 
this section. 
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Figure 2. Monthly and cumulative distributions of observational 
epochs in the 19hr field as of 27th May 2010. There are seasons 
when observations were frequent with a handful of images taken 
every night and there are big gaps when the survey was not sched- 
uled for observation. 



3.1 2D pipeline 

All images taken with WFCAM are processed using an 
image reduction pipeline operated by the Cambridge As- 
tronomical Survey Unit (CASUJE The WFCAM pipehne 
loosely evolved from strategies developed for optical pro- 
cessing (e. g. the Wide Field Su rvey on the Isaac Newton 
Telescope, ijlrwin fc Lewij 20011)) an d implements methods 
originally presented in llrwini ( 19851 ) . The 2D processing is 
a fairly standardized procedure for the majority of projects 
using WFCAM on UKIRT. We give a brief overview of the 
image processing steps here. 

Images are converted into niulti extension FITS for- 
mat, containing the data of the 4 detectors in one pawprint 
as extensions. Other data products from the pipeline (cat- 
alogues, lightcurves) are also stored in binary FITS table 
files. A series of instrumental correction steps is performed, 
accounting for: nonlinearities, reset-anomalies, dark current, 
flatfielding (pixel-to-pixel), defringing and sky subtraction. 
The sky subtraction removes any spatial variation in the 
sky background but preserves its mean level. The sky back- 
ground is calculated as a robust ka clipped median for each 
bin in a coarse grid of 64x64 pixels|j The sky background 
map is filtered by 2D bilinear and median filters to avoid 
the sky level shifting in bins dominated by bright objects. 
The final step is to stack the 9 individual WTS exposures to 
produce one 90 second exposure. Note that a simplified ver- 
sion of the catalogue generation and astrometric calibration 
steps (described below) are run on the individual unstacked 
exposures to ensure that they can be aligned before combin- 
ing. 



^ http : //casu. ast . cam. ac . uk/surveys-proj ects/wf cam 
* The median absolute deviation (MAD) is used as a robust es- 
timator of the root mean square (RMS) in most pipeline compo- 
nents both during processing individual frames and lightcurves. 
For normal distribution, RMS = 1.48 ■ MAD. 



3.2 Catalogues, Astrometry and Photometry 

Object detection, astrometry, photometry and classifica- 
tion are performed f or each frame . Object detection follows 
methods outlined in llrwini (| 19851 ) (see also [Lawrence et al.l 
(2007)). Background-subtracted object fiuxes are measured 
within a series of soft-edged apertures (i.e. pro-rata division 
of counts at pixels divided by the aperture edge). In the se- 
quence of apertures, the area is doubled in each step. The 
scale size for these apertures is selected by defining a scale 
radius fixed at 1.0 arcsec for WFCAM. A 1.0 arcsec radius 
is equivalent to 2.5 pixels for WTS non-interleaved data. 
In 1 arcsec seeing an rcore-radius aperture contains roughly 
2/3 of the total fiux of stellar images. Morphological object 
classification and derived aperture corrections are based on 
analysis of th e curve of growth of the object flux in the series 
of apertures (jlrwin et al.ll2004h . 

Astrometric and photometric calibrations are based 
on matching a set of catalogued objects with the 2MASS 
(jSkrutskie et al.l 120061 ) point source catalogue for every 
stacked pawprint. The astrometry of data frames are de- 
scribed by a cubic radial distortion factor (zenithal poly- 
nomial transformation) and a six coefficient linear trans- 
formation allowing for scale, rotation, shear and coordi- 
nate offset corrections. Head er keywords in FITS files fo l- 
low the system presented in iGreisen fc Calabrettal (|2002h ; 
ICalabretta fc GreisenI (|2002l ). The photometric cal i bratio n 
of WFCAM data is described in iHodgkin etHI (|2009h . 
These calibrations result in the addition of keywords to the 
catalogue and image headers, enabling the preservation of 
the data as counts in original pixels and apertures. 



3.3 Master catalogues 

Following the standard 2D image reduction procedures, cat- 
alogue generation and calibration, we have developed our 
own WTS lightcurve generation pipeline. This is largely 
based on previous work for the Monitor project described in 
llrwin et al.l (|2007l ) where more technical details are given. 

As a first step in the lightcurve pipeline, master images 
are created for each pawprint by stacking the 20 best-seeing 
photometric frames. The master images play dual roles in 
our processing: they define the catalogue of objects of the 
survey for each pawprint with fixed coordinates and iden- 
tifications numbers (IDs). Thus the source-IDs will never 
change for the WTS, however their coordinates could if we 
were to include a description of their proper motions. This 
has not yet been done. 

Source detection and flux measurement is performed for 
each master image to create a series of master catalogues, 
and astrometric and photometric calibration recomputed, 
again with respect to 2MASS. These object positions are 
then flxed for the survey. Each source has significantly bet- 
ter signal-to-noise on the master image, and thus better as- 
trometry (from reduced centroiding errors) th an could be 
achiev ed in a single exposure. As discussed in llrwin et al.l 
(J2007|), centroiding errors in the placement of apertures can 
be a significant source of error in aperture photometry, par- 
ticularly for undersampled and/or faint sources. 

This master catalogue is then used as an input list for 
fiux-measurement on all the individual epochs, a technique 
we call list-driven photometry. A WCS transformation is 
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computed between the master catalogue and each individual 
image using the WCS solutions stored in the FITS headers. 
Any residual errors in placing the apertures will typically be 
small systematic mapping errors that affect all stars in the 
same way or vary smoothly across frames. In practice, this 
can be corrected by the normalization procedure (see Section 
13. 4p . The same soft-edged apertures are used (as described 
above), except that the position of the source is no longer a 
free parameter. Thus for each epoch of observation, a series 
of fluxes for the same sources is measured. 



3.4 Lightcurve Construction and Normalization 

Although the photometry of each frame is calib rated individ- 
ually to 2MASS sources (JHodgkin et al.ll2009l ). these values 
can be refined for better photometric accuracy. Lightcurves 
constructed using the default calibrations typically have 
an RAdS at the few percent level (for bright unsaturated 
stars). To improve upon this, an iterative normalization al- 
gorithm is used to correct for median magnitude offsets be- 
tween frames, but also allowing for a smooth spatial varia- 
tion in those offsets. In each iteration, lightcurves are con- 
structed for all stellar objects and a set of bright stars se- 
lected (13 < J < 17) excluding the most variable decile of 
the group based on their (actual iteration) lightcurve RMS. 
Then for each frame, a polynomial fit is performed on the 
magnitude differences between the given frame magnitudes 
and the corresponding median (lightcurve) magnitude for 
the selected objects. The polynomial order is kept at (con- 
stant) until the last iteration. In the last iteration, a second 
order, 2D polynomial is fitted as a function of image coor- 
dinates. For each frame, the best fit polynomial magnitude 
correction is applied for all objects and the loop starts again 
until there is no further improvement. Multiple iterations 
of the constant correction step help to separate inherently 
variable objects from non-variable ones initially hidden by 
lightcurve scatter caused by outlier frames. The smooth spa- 
tial component during the last iteration accounts for the ef- 
fects of differential extinction, as well as possible residuals 
from variation in the point spread function (PSF) across the 
field of view. 

It is also found that lightcurve variations correlate with 
seeing. In an additional post-processing step, for each ob- 
ject, a second order polynomial is fitted to differences from 
median magnitude as a function of measured seeing. Magni- 
tude values are then corrected by this function on a per-star 
basis. 



3.5 Bad epoch filtering 

Data are taken in a wide range of observing conditions, 
sometimes with bad seeing (~ 2 arcseconds or worse) or sig- 
nificant cloud cover. Additionally some frames are affected 
by loss and recovery of guiding or tip-tilt correction dur- 
ing the exposure. We identify and remove bad observational 
epochs that add outlier data points for a significant num- 
ber of objects in any chip of a pawprint. Where a single 
(corrected) epoch has in excess of 30 per cent of objects de- 
viating by more than 3a from the median flux we flag and 
remove this epoch from all lightcurves. Fig. |3] shows some of 
the per-epoch parameters (which we store in the lightcurve 



files) as a function of 3a outlier object ratio. The number of 
outlying objects has a strong correlation with the residual 
RMS of the second order polynomial normalization (panel a) 
which is a measure of the photometric unevenness of the im- 
age. The (cumulative) frame magnitude offset applied dur- 
ing the normalization (panel c) shows two distinct branches, 
and a large scatter can be seen in the average frame ellip- 
ticity (panel b). These two panels help to identify the main 
causes of bad epochs. High ellipticities arise from frames 
with tracking/slewing problems, while high magnitude cor- 
rections are suggestive of thick (and probably patchy) cloud. 
Our rejection threshold is a compromise between the num- 
ber of affected frames and frame quality. At the outlier ratio 
threshold of 0.3, 39 frames (out of 950, or 4 per cent) are 
removed in the 19hr field. 



3.6 Lightcurve quality 

In Fig. |4] and Table [2] we summarize the WTS lightcurve 
quality at different pipeline steps. In each step we change 
or add one pipeline feature. The normalization procedure 
of the photometric scale by per- frame constant offsets (a), 
the quadratic spatial correction (b) during the normaliza- 
tion, the bad epoch filtering (c) and the seeing correction 
(d) give improvements at the several mmag level for unsat- 
urated bright stars. In panel (d) a theoretical noise model 
curve consisting of Poisson noise, sky noise and a constant 
residual are drawn. A constant systematic error of 3 mmag 
is applied to the model shown in Fig. |3] to bring the model 
roughly into agreement with the data. This should be seen 
as the minimum systematic error in our lightcurves. Some 
saturation appears and makes lightcurve RMS worse for ob- 
jects brighter than J = 13 while for the faint end, the sky 
noise dominates. 



3.7 Transit detection 

To detect transit signals, we use a variant of the boxcar- 
fitting (Box Leas t Squa res, or BLS) algorithm developed by 
lAigrain fc IrwinI (|200J). They derive the transit fitting al- 
gorithm starting from a maximum likelihood approach of 
fitting generalized periodic step functions to the lightcurve. 
It was demonstrated that for planetary transits a simple 
box shaped function is sufficient. The algorithm in this form 
is equ ivalent to BLS developed by lKovacs. Zucker fc Mazehl 
(|2002l ). The signal to red noisqj statistic is used to mea- 
sure transit fitting significa nce (eq. 4 reproduced from 
iPont. Zucker fc Queloj (|2006l )): 



•-^red — 



'Z2i^j ^-i-j 



(1) 



where n is the number of in-transit data points in the whole 
lightcurve, ao is the measurement error of the individual 
data points, d is the fitted depth of the BLS algorithm, dj 
is the covariance of two in-transit data points. Also following 



■^ Noise that includes both uncorrelated (white) and correlated 
(red) components. In some cases, this is called 'pink' noise in the 
literature. 
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Figure 3. Per frame residual normalization RMS (a), average stellar ellipticity (b) and overall magnitude scale correction (c) as a 
function of 3it outlier object ratio in the 19a pawprint. Epochs above a threshold of 0.3 are removed from the survey's lightcurve release 
and candidate search. See text for more details. 



J = 


11-12 


12-13 


13-14 


14-15 


15-16 


16-17 


(a) constant normalization 


9.3 


6.1 


5.7 


6.6 


9.7 


19.1 


(b) quadratic normalization 


9.1 


5.8 


5.0 


5.9 


9.2 


18.8 


(c) + outlier filtering 


8.7 


5.6 


4.9 


5.7 


8.8 


18.1 


(d) + seeing correction 


7.3 


4.8 


4.6 


5.6 


8.5 


17.6 



Table 2. Median lightcurve RMS (mmag) as a function of object magnitudes using different pipeline options for the 19a pawprint. See 
text for details, and Fig. |4]for illustration. 



IPont. Zucker fc Queloa (|2006f ) we adopt a detection thresh- 
old of Srcd = 6. Objects that pass this threshold are consid- 
ered candidate transiting systems. 



4 A SAMPLE OF M DWARFS 

4.1 Identification of M dwarfs 

The high numbers of objects in the WTS and their faint 
magnitudes make spectral classification of all WTS objects 
potentially resource consuming. Instead, we use homoge- 
neous broadband optical (from SDSS DR7. lAbazaiian et al.l 
120091 ')') and near- infrared ("WFCAM) photometry to esti- 
mate reliable stellar spectral types. Specifically, effective 
temperatures are measured from fit ting NextGen stellar 
evolution models (jBaraffe et al.lll99g ) to griz (SDSS) and 
ZYJHK CWFCA M) magnitudes. For comparison, we also 
fit the Dartmouth (jDotter et al.l 120081 ) models, but in this 
case we use 7 passbands (Z and Y are not available). Least 
squares minimization is performed for the available pho- 
tometry, fitting for temperature and a constant magnitude 
offset (distance modulus) as model parameters. The model 
grid data is smoothed by a cubic interpolation. We fol- 
low the temperatu r e-spec tral-class relation in Table 1 of 
iBaraffe fc Chabriei) (| 19961 ): i.e. 3800K, 3400K, 2960K and 
1800K corresponding to spectral types of MO, M2, M4 and 
M9 respectively. 

Fig. E] shows WFCAM and SDSS colour-magnitude 
and colour-colour panels of the stellar objects in the 19hr 
field. The M dwarf sample identified by the SED fitting on 
NextGen magnitudes is marked by red (13 < J < 17) and 
blue (11 < J < 13) crosses. A small number (6 per cent) 
of objects have outlying colour values due to saturation in 



one or more WFCAM or SDSS filters (typically J < 13). 
These magnitudes are sigma-clipped during the SED fitting 
procedure. The panels also show 1 Gyr isochrones from the 
NextGen and Dartmouth models in the 2000K-6500K and 
3200K-7700K temperature intervals respectively. The solid 
and dashed curves correspond to solar metallicity, the dash- 
dot to a metal-rich ([Fe/H]=+0.5), the dash-double-dot to a 
metal-poor ([Fe/H]=-0.5) model isochrone of the Dartmouth 
model respectively. Temperatures of 3800K and 3900K are 
marked (x) on the isochrones. Comparison between Dart- 
mouth model isochrones show little significant difference in 
model colours between ages of 250Myr and 5Gyr. For very 
young stars, we might expect to pick up significant colour 
variation, however we expect very few ve ry young (age < 100 
Myr) stars in our survey field. In fact ICiardi et al.l ([20111) 
analyse the very nearby Kepler field and find that the low- 
mass dwarf population is dominated by young thin disk 
stars, thus our selected age of IGyr is reasonable. 

The figure also presents measured colours of K and M 
dwarfs and Mill giants of th e Pickles photometr ic standards 
(dwarfs: D, giants: v) from ICovev et al.l (|2007l ) and of the 
Bruz ual-Persson-G u nn-St ryker atlas (dwarfs:o, giants:A) 
from iHewett et al.l (|2006l ). All panels show colours in the 
Vega system, AB-y ega offsets are taken from Table 7 in 
iHewett et all (120061 ) . 2MASS - WFCAM conversion s are cal- 
culated by relations given bv lHodgkin et al.l (|2009r ). We de- 
noted by filled green markers the MO dwarf members (one 
and two objects, respectively) of these observations. Our 
identified M dwarfs are separated well from the sample of M 
giants (Vi'^) in the J-H vs. H-K and g-r vs. r-i panels, so we 
expect a low giant contamination level in our sample. The 
panels in Fig.[5]also demonstrate some difficulties in identify- 
ing M dwarfs. Model predictions do not reproduce observed 
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Figure 4. RMS of stellar objects in the 19a pawprint with different pipeline optimizations; a) constant normalization b) quadratic 
normalization c) outlier frame filtering d) seeing correction. In panel d) a noise model (thick red solid line) consisting of Poisson noise 
(dashed line), sky noise (dash-dotted line) and systematic noise of 3 mmag (thin blue solid line) is drawn. 



patterns in all colour combinations. The NextGen isochrone 
has the best agreement in the infrared (J-H vs. H-K) while 
the Dartmouth models fit the r-i vs. i-z colours rather bet- 
ter. We note that the Nextgen colour predictions are too 
blue in the optical bands for low temperatur es, which is a 
know n model attribute (e.g. by 0.5 mag in V-I, iBaraffe et al.l 
Il998l ). 

We also note that based on the residual x^ values during 
the SED fitting procedure the photometric errors from the 
catalogues (shown in the colour-colour panels of Fig[5]) were 
found to be underestimated. We assumed a systematic error 
between the WFCAM and SDSS catalogues and added a 
0.03 mag systematic term in quadrature to the individual 
magnitude errors (used as weights in the fitting). 

The Dartmouth data sets cover a narrower tempera- 
ture range with a lowest temperature of 3200K. This makes 
the dataset unsuitable for selecting later (M4 and later) M 
dwarfs, although it gives a useful comparison for warmer 
stars. We conclude that for early M dwarfs (MO-2), the Dart- 
mouth and Nextgen models give rise to very similar selec- 
tions assuming solar metallicity. In the metal poor and metal 



rich cases the Dartmouth isochrones give about 40 per cent 
increase and 30 per cent decrease in numbers of early M 
dwarfs respectively. 

In the g-r vs. r-i diagram, we note that the scatter in 
g-r for objects with r-i> 1 is larger (by about 0.2 mags at 
g=20) than can be explained solely from the SDSS photo- 
metric errors. This enhanced scatter can be explained by 
reddening alone (see below), and needs no significant spread 
in metallicity. 

iLeggettI (| 19921 ) analyzed photometry of 322 M dwarfs 
and found that the effects of metallicity can be seen in M 
dwarf infrared colours (I-J, J-K, I-K, J-H, H-K) while not 
discernible in visible colours (U-B, V-I, B-V). They also 
noted tha t this feature i s not reprodu ced by evolutionary 
models of iMouldl (| 19761 ): lAUardl (|l990l ). We note that the 
D artmouth mode l also shows a much smaller effect than seen 
in iLeggetd (| 19921 ) in the infrared, but a very large effect in 
g-r vs r-i. 

In our paper, we base object classification on the 
NextGen model. Out of the 59000 morphologically identi- 
fied stars in the 11 < J < 17 magnitude range, we identified 
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Figure 5. Colour-magnitude and colour-colour plots of the M dwarf sample (13 < J < 17, red crosses; 11 < J < 13, blue crosses) of this 
study and other stellar objects (10 < J < 16) in the 19hr field bas e d on WFCAM and SDSS data. A sequence of spectroscopic K and M 
dwarfs are shown from lCovev et alj 1120071 ) (D) and lHewett et al.l ll2006f) (o), respectively. Green markers show the MO dwarf members 
of these sequences. Triangles (v,'^) are M giants (III) from the same studies. Estimated maximum interstellar extinctions in different 
colours are marked by red arrows. Error crosses are from the corresponding catalogues. 



~ 4600 M dwarfs. A more in-depth modelling of the objects 
in the WTS is an ongoing effort. 



4.2 Interstellar extinction 

During the SED fitting procedure, we cannot fit for the 
interstellar reddening and for the distance modulus on a 
per object basis due to being the fit badly determined. 
Therefore, we make an upper estimation for the redden- 
ing in our sample. We estimate a maximum distance for 
target M dwarfs of 1.5 kpc by assuming that one of 
our intrinsically brightest objects (an MO, Mj — 6 from 
NextGen) is observed at the faint limit of the study (J = 
17). Ay ext inction values are calculated from the Galac- 
tic mo del of iDrimmel. Cabrera-Lavers fc Lopez- Corredoiral 
(12003). The three dimensional Galactic model consists of 
a dust disk, spiral arms mapped by HIT regions and a lo- 
cal Orion-Cygnus arm segment. They use COBE/DIRBE 
far infrared observations to constrain the dust parame- 



ters in the model. We use the code provided to deter- 
mine Av at 1.5 kpc. To calculate absorption in UKIRT 
and SDSS bandpasses, we use conversion f actors A/ Ay from 
Table 6 of ISchlegel. Finkbeiner fc David (Il998l) that evalu- 
ates the reddening l aw oflCardelli. Clayton fc Mathid (| 19891 ) 
and of lO'Donnelll (|l994l ) for the infrared and the opti- 
cal bands respectively. The reddening laws assume Rv = 
Av/Eb-v ~ 3 .1, an average value for the diff use inter- 
stellar medium (|Cardelli. Clavton fc Mathid 1 19891 ). Red ar- 
rows show our reddening estimates in different colours in 
Fig.\E\{Av : 0.409, Aj: 0.113, Ej-k: 0.067, Ej-h- 0.041, 
Eh-k: 0.026, Eg-r: 0.130, Er-i: 0.083, Ei-,: 0.076). The 
total Galactic extinction in t he 19 hr field direction from 
ISchlegel. Finkbeiner fc David (|l998D is Ay = 0.43^ 



^ Provided by the Nasa 
http : //ned . ipac . caltech . edu/ 



Extragalactic Database. 
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4.3 M dwarf bins 

Stellar radius changes significantly between early- and late- 
type M dwarfs. For our sensitivity simulation purposes, we 
can use a linear approximation for the mass-radius relation- 
ship for M dwarfs: M^/Mq ~ R,/Rq. This sli ghtly deviates 
from the NextGen mass-radius prediction (jBaraffe et al.1 
Il99a ) b ut it is in good ag reement with observed M dwarf 
radii in iKraus et al.l (|201ll V It is important in our simula- 
tions to treat early M dwarfs distinctly from later smaller 
M dwarfs. Therefore we divided our M dwarf sample into 
three coarse spectral bins. In the 19hr field, 2844 stars were 
identified as MO-2 (3800K > Teff > 3400K), 1679 as M2-4 
(3400K > Tcft > 2960K) and 104 as M4-9 (2960K > T^s > 
1800K). The third, latest (M4-9) type bin is very sparsely 
populated (2 per cent of the M dwarf sample) and numbers 
suffer from high uncertainty, so this bin will be omitted from 
the simulations, leaving us with two subsamples. 

Considering the 'colour distance' along the isochrones 
between the boundary of our identified object groups, the 
Pickles MO members (green squares) and the dependency 
of the isochrones on model parameters, we estimate the bin 
temperature edges to be uncertain about 250K. 



4.4 M Dwarf Lightcurves 

M dwarfs are known to be intrinsically more variable 
than more massive ma in sequence stars prim arily due to 
spots (e.g. IChabrier. Ga llardo & Baraffc 2007,: ICiardi et al.l 
l2011f ). The fraction of active M dwarfs also rises towards 
later (M4-9) subtypes (|West et al.l 120111 ). We examine the 
variability of the M dwarfs compared to the rest of the stel- 
lar sample in Table |3] and Fig. [S] The photometric precision 
is at the 4-5 mmag level for brighter objects (J < 15) and 
at the percent level for the fainter region (15 < J < 17, 
86 per cent of the M dwarfs). We compare the binned me- 
dian RMS values of the three M subtypes to those of the 
warm (T>3800K) stars. Although the effect is small, we do 
find evidence that the median RMS for all three M dwarf 
samples is systematically higher, although only at about the 
0.5-1 mmag level, at all magnitudes (see Fig|5J. 

We also calculate the fraction of 3a outliers (a proxy for 
the most variable objects) for the different subsamples. We 
find that all three M dwarf samples, as well as the warmer 
stars exhibit a ten per cent outlier fraction. 

Finally, these findings are also supported by 2D 
Kolmogorov-Smirnov tests (JFasano fc Franceschinil [l987J). 
The RMS versus magnitude distributions are found to be 
consistently different for each of the M dwarf samples ver- 
sus the warm stars to very high confidence (test p values 
< 10~^). These results add support to the argument that 
observing in the J b and gives limited sensitivity to spots 
l|Goulding et al.ll2012l '). 

Typical signal depths for edge-on planetary systems are 
also shown in Fig. [6] to give a rough impression of our sensi- 
tivity. We expect to be able to detect Jupiter size planets in 
edge-on systems around all of our M dwarfs (in all spectral 
bins: MO-2, M2-4, M4-9), and Neptunes only in the M2-4 
and M4-9 bins. As stated above, the latest type bin (M4- 
9) is sparsely populated, and omitted from the remaining 
analysis. 



J = 


12-13 


13-14 


14-15 


15-16 


16-17 


MO-2 


4.4 


5.1 


5.7 


8.5 


17.4 


M2-4 


4.8 


4.8 


6.4 


8.5 


17.5 


M4-9 


3.9 


5.6 


4.8 


9.9 


17.8 


Earlier types 


4.8 


4.5 


5.4 


8.3 


16.6 



Table 3. Median RMS (mmag) of M dwarf lightcurves and of 
earlier dwarfs for different object magnitudes. 



4.5 The WTS sample compared to other surveys 

The WTS is not the only transit survey specifically aimed 
at discovering planets around M dwarfs, and it is useful to 
place our survey in context. The MEarth survey individually 
target s about 2 000 bright M dwarfs in a custom i + z-hand 
filter (JBerta et al. 2012 ). The survey was designed to focus 
on the nearest and brightest mid M dwarfs and has dis- 
covered one Super Earth size planet at the time of writing 
(|Charbonneau et aLll2009h . 

The Palomar Transient Factory (PTF, iLaw et ali 

(|2012l )) is a project targeting M dwarfs in the R band. Their 
goal is to observe a total of 100,000 M dwarfs during the sur- 
vey lifetime collecting about 300 epochs in each observing 
period. The PTF operates in roughly the same magnitude 
range (specifically over R=14-20) as the WTS and has a 
similar, seasonal observing schedule. 

The Pan-Planets project also focuses on late-type stars. 
Their goal is to target approximately 100,000 M dwarfs 
in 40 square degrees. The programme collected about 100 
hours of observations using the Panoramic Survey Tele- 
scope and Rapid Respon se System (PanSTARRS) so far 
JKoppenhoefer et al.llJOOgl ). 

The Next- Generation Transit Survey (NGTS) is also 
an upcoming initiative to target a high number of bright 
M dwarfs (y<13) using wide field cameras. The prototype 
instrument observed ~ 400 early and ~ 50 late M dwarfs 
in 8 square de grees FOV at the on e per cent photometric 
precision level (jChazelas et al.ll2012h . 

Regarding space missions, the Kepler Mission in its 
Q2 data release has a relatively small (1086) number of M 
dwarfs but their data quality is much higher than for ground 
based programmes. The mission also samples uniformly in 
time with no large data gaps. We briefly discuss the planet 
occurrence rate in this sample in Section [7] based on IH12| 
and compare it with our results. 

Compared to these other transit surveys, the WTS is 
targeting a reasonable sample, even using only one of the 
four fields (at the time of writing three of the four fields 
are lagging in coverage). The PTF and Pan-Planets samples 
are potentially ground-breaking, if they can obtain enough 
observations. 



4.6 A lack of hot Jupiters discovered around M 
dwarfs in the WTS 

In the WTS, the initial candidate selection criterion 
is currently a well-defined decision based solely on the 
signal-to-noise ratio. Exhaustive eyeballing and photomet- 
ric/spectroscopic follow-up of the small numbe r of candidate 
planets around M dwarfs (|Sip6cz et al.l 120131 ) has revealed 
them all to be grazing eclipsing binary systems, or other 
false-positives. This search was complete for all --^Jupiter- 
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Figure 6. Comparison of liglitcurve RMS in tlie M dwarf spectral groups. The solid curve shows median values of the RMS per magnitude 
bin. The bottom right panel shows the differences of the median RMS model curves from the median RMS of the earlier type stars. {MO-2: 
solid, M2-4; dashed, M4-9: dash dotted). Signal depths of edge-on transiting systems with Neptune, Jupiter or Super-Earth (2.7-R0) size 
planets are marked. 



sized candidates down to J=17 in the 19 hour field, although 
avoiding objects with periods very close to one day. Two hot 
Jupiters h ave however been discovered around mor e mas- 
sive hosts (|Cappetta et al.ll2012l : Isirkby et aLll2013h . This 
confirms that the WTS is sensitive to transits at the ^ 1 
per ce nt level (and indeed t hese two objects are very ob- 
vious). IZendeias et al.l (|2013l ) have independently produced 
lightcurves (using difference imaging) and searched for can- 
didate planets in the same data, using stricter autoinated 
classification, and also find no candidate Jupiters around M 
dwarfs in the 19 hour field. This lack of hot Jupiters discov- 
ered around small stars in the WTS is what motivates the 
rest of this study, to address the question: how significant is 
this result. 



5 SIMULATION METHOD 

5.1 Recovery ratio 

We determine the sensitivity of the WTS for four distinct 
star-planet scenarios. For the two M dwarf subsamples, we 



No. Sp.bin Teff (K) M,IMq R./Rq M^/Mj^p fip/fl.Tup 



1 


M2-4 


2960-3400 


0.4 


0.4 


1.0 


1.0 


2 


M2-4 


2960-3400 


0.4 


0.4 


0.054 


0.36 


3 


MO-2 


3400-3800 


0.6 


0.6 


1.0 


1.0 


4 


MO-2 


3400-3800 


0.6 


0.6 


0.054 


0.36 



Table 4. Simulated planetary systems 



consider Neptune and Jupiter size planets in short period 
(0.8-10 days) orbits around early (MO-2) and later type (M2- 
4) stars. The mass and radius of the star and the radius 
of the planet are kept fixed in each scenario (see Table |4]). 
We use conservative assumptions for the stellar parameters: 
stellar radii are overestimated, using the maximum values 
in each spectral bin (i.e at MO and M2), and planet radii 
are (under)estimated assuming solar system radii even for 
hot planets. The expected number of recovered planetary 
systems can be written as: 



Ni 



-/VstarsJ-Pdct 



(2) 
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where A^stars denotes the number of stars in the actual star 
subtype group, / is the (unknown) fraction of the stars that 
harbour a planetary system, Pdct is the (average) probabil- 
ity of discovering a system by the survey around one of its 
targets if we assume that the star harbours a (not neces- 
sarily transiting) planetary system. Pdot is expressed as a 
function of planetary r adius (Rj,) and orbit al period (T). In 
general, we can write (jHartman et al.ll2009l ): 



-Pdct — 



I I PrPTp{Rp,T)dRpdT 



(3) 



where Pr is the average recovery ratio, i.e. the average (con- 
ditional) probability of recovering a transit from a lightcurve 
if the lightcurve belongs to a transiting system, Pt is the 
geometric probability of having a transit in a randomly ori- 
ented planetary system, and p{Rp, T) is the joint probability 
density function of Rp and T for planetary systems. We de- 
ter mine these terms in E g. [3] separately in each scenario. See 
also lBurkeet all (|2006l '): lHartman et al.l (|2009l 'l. 

We consider circular orbits only (e=0). Pt can be 
given analytically, for a random system orientation, Pt = 
(-R* +Rp)/a where a denotes the semi-major axis of the sys- 
tem, R, is the stellar radius. The joint probability density, 
p{Rp,T), must be treated as a prior. We consider planetary 
configurations at discrete Rp values only, so the Rp depen- 
dence of the density function simplifies to a (5-function. As 
the period dependence is ill-constrained, we will discuss uni- 
form and power-law functions as a prior distributions. In the 
IH12| study, a power law model with a n exp onential cutoff at 
short periods was fitted (Table 5 in lH12r ) to Kepler giant 
planet detections around mostly GK dwarfs (see sample cri- 
teria in Section [7|. We use their model function normalized 
to our studied period range as prior distribution. 

Pr is determined numerically by the Monte Carlo iter- 
ation loop. 



5.2 Determining Pr 

We identify a set of 4700 quiet lightcurves in the 19hr 
field that serve as input for the simulations. This simula- 
tion sample consists of mostly M dwarfs (and slightly hot- 
ter stars) covering the magnitude range 11-17. The sam- 
ple preserves the observed distribution of M dwarf apparent 
magnitudes in the WTS. By quiet we mean that the un- 
perturbed lightcurves show no signature of a transit (BLS 
Srcd <= 6). By adding noise free signals to these lightcurves 
and recovering them from the noisy data, we can quan- 
tify the effect of noise on Pr. As shown in Section [4.41 we 
found little difference between the noise properties of the 
lightcurves for the MO-2 and M2-4 subclasses. We simulate 
large numbers of transiting exoplanet systems to determine 
the recovery ratios for the four scenarios under investigation. 
In each iteration, a transiting planetary system is created 
with parameters randomly drawn from fixed prior distribu- 
tions as detailed below (assuming a circular orbit). A simu- 
lated transit signal is then added to the randomly selected 
(real) lightcurve. We try to recover the artificial system us- 
ing the transit de t ection algorithm discussed in Section [3.71 
(jAigrain fc IrwinI (1200J)). Pr is estimated as the ratio be- 
tween successful transit recoveries and the total number of it- 
erations. We make a distinction between two different cases. 
In the threshold case we consider the signal successfully re- 



covered if the detection passes the same signal to red noise 
level as used for WTS candidate selection (S'rcd>6). In the 
periodmatch case, we additionally require that the recovered 
period value matches the simulated one. This is discussed in 
more detail in Sec. 16.31 

The period value (T) is drawn from a uniform distribu- 
tion in the range of 0.8 to 10 days. The period determines 
the semi-major axis (a) of the system, given the masses 
of the star (assumed to be 0.6 or O.4M0) and the planet 
(assumed to be 1.0 or 0.054 Mjup). A randomly oriented 
system is uniformly distributed in cos i where i is its or- 
bital inclination. The random inclination is chosen to satisfy 
^ cos j < (_R* -I- Rp)/a to yield a transiting system. This 
also allows for grazing orientations. The phase of the transit 
is also randomly chosen from within the orbital period. 

Observed dates in the target lightcurve are now com- 
pared with predicted transit events. If there would be no 
affected observational epochs, then the iteration ends and 
the generated parameters recorded. Otherwise, a realis- 
tic, quadra tic limb dark ening model is used with coeffi- 
cients from IClaretl ([20001) to calculate brightness decrease 
at in- transit observational epochs (Mandcl & Agol 2002; Pall 
l200g ). This artificial signal is added to the lightcurve magni- 
tude values, and BLS is run on the modified lightcurve. Both 
generated and detected transit parameters are recorded for 
the iteration. We note that the transit detection algorithm 
is the most computationally intensive step in the loop. A 
total of 75,000 iterations were performed. 



6 RESULTS AND DISCUSSION 

6.1 Sensitivity effect of the observation strategy 

The fiexible observing mode of WTS has an inherent limi- 
tation on our sensitivity to short-period transiting systems, 
and it takes multiple seasons to build up enough epochs to 
reliably detect them. In Fig.[7]a simple sensitivity diagram is 
shown which considers only the actual distribution of obser- 
vational epochs for the 19 hour field. We use the simulated 
transiting systems from the Neptune-size planets around an 
MO dwarf scenario and compare the simulated transit times 
with our real observational epochs. A system is considered 
detectable here if at least 5, 10 or 15 in-transit observational 
epochs occur calculated from the simulated parameters. The 
fraction of detectable systems has an obvious strong depen- 
dence on the period of the transiting system and the re- 
quired number of in-transit observational epochs as well. Of 
course, it depends on the noise properties of our data how 
many in-transit observations are necessary for detecting a 
transit event. 



6.2 Detection statistic 

We show the (signal to red-noise) detection statistic distri- 
butions in Fig. (5] The black curve belongs to the survey M 
dwarfs (unmodified data) in the 19hr field (selection criteria 
described in l4.3|l . The coloured curves are for the simulated 
transiting systems for the four scenarios. For comparison, 
they are normalized to the total number of M dwarfs in the 
survey, i.e as if every WTS M dwarf target were either an MO 
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Figure 7. The effect of the observation strategy on the WTS 
sensitivity to short period transiting systems in the 19 hour field. 
We use the actual epochs of our 19 hour field observations com- 
bined with a large sample of simulated planets. The fraction of 
simulated transiting systems is shown where at least 5, 10 or 15 
actual WTS observational epochs coincide with the simulated in- 
transit times. This illustration does not use the noise-properties 
of the WTS, thus does not consider the case as to whether the 
events are actually detected or not. 



or an M2 with a transiting Jupiter or a Neptune (which may, 
or may not, have transited during the WTS observations). 

The original (unperturbed) WTS M dwarf S'rcd distribu- 
tions shows a marked tail above the detection threshold (25 
per cent of all objects, see Sec l6.4[l . This could be caused by 
effects such as: correlated noise, intrinsic variability (spots), 
eclipsing binaries or possibly transiting planets. The simula- 
tions show that for initially "quiet" lightcurves (S'rcd < 6) the 
injection of a planet signal (transiting, but could be graz- 
ing) can perturb the measured signal-to-noise value above 
the detection threshold. (M2+J:50, M0+J:39, M2+N: 5.1, 
MO+N: 4.1 per cent of all iterations) In other words, the 
WTS is sensitive to Jupiters, and rather less sensitive to 
Neptunes. Given the other possible causes for a perturbed 
Srcd, the shape of this distribution does not in itself present 
a direct measurement of the planet population. It is worth 
noting that for the Jupiters (M0+J,M2+J), the recovered 
Sred values can be rather higher than the largest values we 
actually see in the WTS. 

The Neptune cases have much smaller residual tails 
which sit only a little above the detection threshold. These 
signal injections cause detection statistic increases above the 
detection threshold only in a small number of cases. There is 
also little dependence on the size of the host star (the green 
and red curves are similar). Our simulations, with the same 
stellar magnitude distribution as the survey, are dominated 
by fainter objects. The similarity of the detection statistic 
distributions in the figure indicates that for Neptune sized 
planets, the survey sensitivity depends on the stellar radius 
only for the brightest stars. For fainter objects, the sensitiv- 
ities are the same and as we see later, they are equally low 
(Figures [TOl and [TT|) . 



6.3 Detected period 

While passing the signal-to-noise statistic threshold in the 
simulation is caused by the injected signal, the detected (re- 
covered) parameters, particularly the period, may not be 
correct. Indeed, it is common in ground based surveys to find 
multiple peaks in the BLS periodogram containing aliases of 
the observing window function. 

The detected period is perhaps the most important of 
the model parameters: At the eyeballing stage of candidates, 
this period has the most influence in judging a transit detec- 
tion real or false. In a lightcurve folded on the correct period 
(or a harmonic thereof), the transit signals are visually eas- 
ily recognized. When folded on a random period, in-transit 
points are hardly distinguishable from random outliers. A 
good initial period value is also important for the timing 
of follow-up observations. With this in mind, we try to ac- 
count for the importance of this effect in our simulations. 
We consider two recovery rates: 

• threshold — Srcd exceeds 6.0, and 

• penodmatch — we additionally require that the de- 
tected period value matches (or be a harmonic of) the gen- 
erated period. We allow factors of 1, 5/4, 4/3, 3/2, 5/3, 2, 
5/2, 3, 4, 5 between the two values. 

This step uses external information and obviously we 
cannot know the true period in the actual survey. However, 
in the simulations, it does allow us to place more strict cri- 
teria on detected simulated planets, and mimic the effects 
of differentiating between high priority candidates, which 
are likely to receive follow-up time, and low-priority candi- 
dates for which follow-up may be too expensive. This effect 
is much more significant for the Neptunes, where individual 
events are shallow. 

We note that the number of in-transit data points is 
correlated with recovering the correct period. In the M2+J 
case, we find that if a simulated system has 10 in-transit data 
points we detect the correct period 80 per cent of the time, 
reaching almost 100 per cent with 20 in-transit data points. 
For larger stars (shallower transits), more in-transit data 
points are needed (MO+J: 12 and 30 points for 80 and 100 
per cent respectively). For Neptunes even more transits are 
required to secure the period. For the M2+N case, we find 20 
in-transit points recovers 50 per cent of the detected systems 
with the correct period. For the MO+N case the situation is 
even worse, and we almost always detect an alias (see Fig.|5}. 
This confirms that the BLS algorithm is "lost" at these low 
signal levels. Comparing the necessary number of in-transit 
data points for reliable transit detections, according to Fig. 
[T] we can see an inherent limitation of the transit searching 
efficiency of the WTS and possibly other ground based low 
cadence surveys (e.g. PTF). 

In this study we determine Pr for the threshold and pe- 
riodmatch cases. The real value of Pr that characterizes the 
WTS probably lies between these two extremes and depends 
on all the (sometimes subjective) steps of the candidate se- 
lection and follow-up strategy. Having fewer quality criteria 
in the survey increase the number of candidates at the cost 
of higher false alarm ratio and more follow-up resource con- 
sumption. 

We note that surveys with close to real time data pro- 
cessing and instant access to follow-up observation facilities 
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Figure 8. Detected signal to red noise histograms (bin size is 0.2) of the transit injected quiet lightcurve iterations (coloured) for Jupiter 
(left) and Neptune (right) size planets in the four simulated scenarios. For comparison, the unmodified (black) M dwarf lightcurves of 
the survey is shown and curves of the iteration loops are normalized to the number of M dwarf objects in the survey (black curve). The 
vertical line indicates our candidate selection threshold of 6. 



could use a different approach and need not necessarily bal- 
ance between quality of recovered transits and false positive 
rates. They can maximize the probability of detecting ac- 
tual transit events by selecting follow-up times optimized 
to their current data and iteratively revise predictions as 
data accumula tes. A Bayesian approac h to such a strategy 
is described in lDzigan fc Zucken (|201ll ). 

In Fig. [9] scatter diagrams and histograms of simu- 
lated and detected periods are shown for recovered simu- 
lated planets. We compare our most sensitive case in terms 
of injected signal depth (M2+J), to the least sensitive case 
(MO+N). They show a clear contrast in the quality of de- 
tected (signal to noise selected) transit signals. 

For Jupiters, about 70 per cent of the detections also re- 
cover the injected period or its harmonic value (middle row 
in Fig. [9]). Strong harmonic lines can be seen in the scatter 
panel (top left). For MO stars (not shown), the ratio of recov- 
ered period values is even better. The shallower transits are 
longer in duration for a given period. The detected periods 
(bottom row, filled histogram) are biased towards shorter 
values peaking around 3 days, with aliases appearing at 1 
and 2 days. In these iterations the injected signal changes 
the lightcurve in a way that aliased periods give the most 
significant box-fit signals. We see a gently decreasing trend 
for the periods of simulated systems that were recovered 
(bottom row, empty histogram), i.e. we are less sensitive 
to longer period planets, and more likely to underestimate 
their periods. We can conclude that signals of Jupiter sized 
planets can firmly be detected around all of the M dwarfs 
in the survey with good initial period values. 

The Neptunes are shown in the right-hand column of 
Fig. [9] In accordance with our qualitative assessment based 
on the lightcurve RMS, signals of Neptune sized planets are 
at the boundary of being lost in the lightcurve noise in the 
WTS. In the M2+N case, only 12 per cent of the detec- 
tions have also matching periods, in the MO+N case the 
correct period is almost never recovered by the box fitting 
algorithm (middle row) . Recovered period values are heavily 



dominated by alias values. We conclude that Neptune sized 
planets can be detected in the survey only in favorable cases 
around smaller (later) M dwarfs. 



6.4 False positives 

In a transiting planet survey, we have to deal with two types 
of false positives at the initial transit candidate discovery 
step. The detection statistic may pass where the transit box 
model is fitted just on (i) random fiuctuations or systemat- 
ics. In this case, the transit signal does not exist at all. (ii) 
Physical signals in the lightcurve may also belong to various 
eclipsing binary configurations or to variable stars (e.g. from 
spots) which can mimic transit events for the algorithm. 

The amount of additional analysis to rule out false posi- 
tive candidates can vary from little additional manual check- 
ing, through refinement of transit parameters, up-to obtain- 
ing additional observational data with better cadence. Graz- 
ing binary stellar systems can mimic planetary transits be- 
yond the survey's photometric precision and require addi- 
tional spectroscopy measurements to be ruled out with high 
confidence. 

The false positive ratio of the candidate selection pro- 
cedure of the survey cannot be quantified by the simula- 
tion alone. A lightcurve passing the detection threshold with 
known injected simulated transit signal cannot be a false 
positive per se. 

So far, all the actually selected, eyeballed and later 
foUowed-up (by additional photometry and/or spectroscopy) 
candidates aro und M dwarf s turri ed out not t o be a plan- 
etary system (|Sipocz et al.l 120131 ). Following iMiller et al.l 
1)20081 ) ■ this can be used to estimate the false positive ra- 
tio of the survey. This encompasses both false positive cases 
above. We assume that all Sicd > 6 detections in our original 
(unmodified) lightcurves are false positives. At 25 per cent 
(black curve above the marker in Fig. [HI), this is a rather 
large fraction. This number gives a simple (worst-case) de- 
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Figure 9. Properties of period values of iterations that pass the signal-to-noise threshold in the most (M2+ Jupiter, left) and least 
(MO+Neptune, right) sensitive scenarios. Best fit detected period values as a function of simulated period values (top row) and normalized 
histograms of simulated/detected period ratios (middle row; bin size is 0.05). Bottom row: Normalized histograms (bin size 0.1) of detected 
(filled) and simulated periods (empty) i.e. projections to axes of top row scatter panels. 
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scription of our false positive ratio, and cannot be usefully 
applied to the simulations. 



6.5 Transit recovery ratios 

We show results of Pr as a function of stellar magnitude for 



0.8 - 3.0 



3.0- 5.0 



5.0 - 10.0 



We adopt 0.8- 



Hartman et al.l 



ranges of simulated period in Figs. llOl andfTT 

3.0, 3 .0-5.0, 5.0-10.0 day ranges following 

((20091) for extremely hot Jupiters (EH J), very hot Jupiters 
(VHJ), and hot Jupiters (HJ) respectively. Filled red circles 
(•) and green crosses (+) represent the threshold and period- 
match cases respectively. Integrated detection probabilities 
(fdet) weighted by different prior assumptions including the 
geometric probabilities of transiting orientations are shown 
in Table El 

For the Jupiter cases (upper rows in both figures), 
the WTS sensitivity has a maximum around J=13.5 and 
drops towards fainter (J>15-16) objects. This is in accor- 
dance with our expectations, we have higher noise levels 
towards fainter objects but there is occasional saturation at 
the bright end (J < 13). There is little dependence on stel- 
lar radius, appearing only for the fainter stars. The threshold 
curves are less affected by simulated period than the peri- 
odmatch recovery rate. In the shortest period window the 
threshold and periodmatch values are practically the same. 
In the longest period panels the threshold ratio is about twice 
that of the periodmatch ratio showing that while the signal- 
to- noise detection statistic can recover signals, many of these 
systems may be missed due to poor initial period guesses. 
The survey has a much lower sensitivity for Neptunes (bot- 
tom rows in Figs. [10] and llip . The threshold curves show 
roughly 25 per cent of the recovery rate seen for Jupiters 
around bright stars. The difference between the threshold 
and periodmatch case is also more significant. The period- 
match recovery rates are close to zero for all MO cases and 
very low for the M2 cases as well. The only exception is the 
best signal-to-noise case around M2 stars, in the short period 
window for bright objects where threshold rates reach half of 
the Jupiter value and the periodmatch rates are not close to 
zero. It looks like the WTS survey only really has a chance 
of discovering extremely hot Neptunes around late (M2 or 
later) and bright M stars (J < 15). The lack of an accurate 
period determination would make follow-up of Neptune size 
candidates difficult in all other cases. 
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Figure 10. Recovery ratios (Pr) of simulated transiting system 
for M0+ Jupiter (upper row) and MO+Neptune (lower row) sce- 
narios in three period ranges as a function of stellar brightness. 
Filled circles (•) and crosses (+) represent the threshold and pe- 
riodmatch ratio respectively. 
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Figure 11. Recovery ratios {Pr) for M2+ Jupiter (upper row) 
and M2+Neptune (lower row) simulated scenarios. 



6.6 Error considerations 

We use the seeing corrected lightcurves as a source for our 
simulation. While it could be more accurate to modify mea- 
sured raw fiux values and to run all the pipeline components 
to create the modified lightcurves, it would be unreasonably 
resource consumin g, nor assumed to have an impact on our 
sensitivity results. iBurke et al.l (|2006l ) decrease the sensitiv- 
ity with a constant estimation of 0.03 to deal with their trend 
filterings' fiattening effect. In our case it is assumed to have 
a negligible difference only as we do not use trend filters that 
uses solely lightcurve data. As it is described in the pipeline 
overview, images are transformed into lightcurves simultane- 
ously, calculating magnitude scale offsets using non-variable, 
bright objects on each frame. Compared to the several thou- 
sand objects per frame, only a small fraction of objects may 



have a transit on a given frame thus the determined magni- 
tude offset value would be independent of an injected signal. 
Being a bad weather fallback project, the survey typically 
has a few (2-3) observational epochs per observing night. 
Thus we do not expect significant correlation between the 
seeing and the timing of transit signals. 

6.7 Limits on planetary fraction 

Assuming that the survey planetary candidates are all false 
positives, as our follow-up suggests (see Sec. I4.6|l . we can 
place upper limits on the planetary fraction (/) for our four 
test cases. The assumption of complete follow-up is much 
less secure for Neptune-sized planets, nevertheless the upper 
limits remain valid for the survey until a signal is detected. 
While Eq. [2] gives the expected number of detections, the 
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System type 


prior 


A^stars 


Pdct 


/95% 


Pdot 


/95% 








sn 


sn 


pm 


pm 


M0+ Jupiter 


Kep. 


2844 


0.0363 


2.9% 


0.0311 


3.4% 


M2+Jupiter 


Kep. 


1679 


0.0414 


4.3% 


0.0346 


5.2% 


MO-4+Jupiter 


Kep. 


4523 


0.0382 


1.7% 


0.0324 


2.0% 


M0+ Jupiter 


uni. 


2844 


0.0360 


2.9% 


0.0305 


3.5% 


M2+Jupiter 


uni. 


1679 


0.0405 


4.4% 


0.0332 


5.4% 


MO-4+Jupiter 


uni. 


4523 


0.0377 


1.8% 


0.0315 


2.1% 


MO+Neptune 


Kep. 


2844 


0.0027 


39% 


0.0001 


100% 


M2+Neptune 


Kep. 


1679 


0.0025 


71% 


0.0003 


100% 



Table 5. Planet detection probability and upper limit results on 
planetary fractions in the WTS in the signal-to-noise threshold 
(sn) and periodmatch (pm) interpretation. The upper part shows 
limits for the simulated Jupiter scenarios (M0,M2) and its inter- 
pretation for the wh ole M O-4 dwarf sample. Values are shown for 
priors based on the IH12| study (Kep.) and uniform (uni.). The 
lower part shows representative values for the Neptune cases. 



by the number of M dwarfs in the corresponding groups. 
The WTS can put an upper limit of /950/, =1.7-2.0 (1.8- 
2.1 for uniform) per cent on the occurrence rate of short 
period Jupiters around early-mid (MO-4) M dwarfs. We dis- 
cuss these results in context in Section [T] 

For Neptunes, the WTS detection probabilities are an 
order-of-magnitude smaller than for Jupiters. We also find a 
strong dependence both on the assumed prior for the period 
distribution and the threshold/ periodmatch interpretations. 
As seen, the recovery of Neptunes in the WTS is a difficult 
and unreliable task. The probability of recovering a Neptune 
with the correct period [periodmatch Pdet) is so low because 
detected signals are dominated by aliased periods. Thus it 
is very difficult to distinguish genuine Neptunes from false 
detections. The main way to improve on the WTS sensitivity 
to Neptunes would be to improve on the noise properties in 
the data. This is beyond the scope of this paper. 



actual number of detections follows a Poisson distribution. 
If the expected number of planets is A^dct , the probability of 
detecting k planets is; 



Pk = 



k\ 



exp(-iVdct) 



(4) 



We use this expression at fc = as a likelihood function for 
A^det- We require that A^dot be within our ^ A''dct < A^max 
confidence interval with 95 per cent confidencqj- Solving 



P{Ni^t < A^max) = 



JVmax 

I ex.p{- N)dN = 0.95 



(5) 



for A^max, we get Ai'max = 3.0. 

A^dct < 3.0, using Eq.[2] we get: 



/95% ^ 



3.0 



A^.i 



From the requirement of 



(6) 



Assuming zero detections and applying the above tech- 
nique, we can compute robust upper limits on the fraction 
of planet host stars in our sample. We show the threshold 
and periodmatch Pdct values along with the planetary frac- 
tions in Table [S] We show results for the two different as- 
sumptions for the period dis tribut ion: best fit exponential 
cut power law from Kepler (|H12| ) (Kep.), uniform (uni.). 
For the hot Jupiters as a whole (periods 0.8-10 days), we 
can place an upper limit of /95% =2.9-3.4 per cent planet 
fraction around the early M subtype group (MO-2) assum- 
ing the Kepler power law prior period distribution. Using 
the simple uniform prior gives similar sensitivities (2.9-3.5 
per cent) in this short period range. The smaller sample of 
cooler stars leads to a higher upper limit of 4.3-5.2 (4.4-5.4 
for uniform) per cent for the M2-4 group. For Jupiter size 
planets, the threshold and periodmatch detection probabil- 
ities are essentially the same, and the chosen prior period 
distribution has also little impact. 

We can combine the M0-M4 stars into one bin by calcu- 
lating the average sensitivity for the whole sample weighted 



^ This is actually a Bayesian interpretation for the 95 per cent 
probability credible interval for the parameter of the Poisson dis- 
tribution, assuming a sufficiently wide, uniform prior for the pa- 
rameter. For nonzero cases in Sec[71 intervals with equal posterior 
probability values at the endpoints are chosen. 



7 DISCUSSION: COMPARISON WITH 
KEPLER AND RV STUDIES 

Although Kepler is not a specialized survey for M dwarfs, 
due to its exceptional quality data and large number of tar- 
gets, it is probably the highest impact exoplanetary transit 
survey to date. In this section we compare our result s to t he 
planet occurrence study of the Keple r data in paper IH12| . 

In the second part of the IH12| study, the planet oc- 
currence rate is determined in three planet radius bins (2- 
4_R0, 4-8i?0, 8-32_R0) as a function of stellar type (effective 
temperature; in 500K wide temperature bins from 3600K to 
7100 K; see figure 8 inlHll) . They use the Q2 Kepler data re- 
lease IjBorucki et al.ll201ll ). In each stellar temperature bin, 
the total planet occurrence rate is calculated by adding the 
contribution of each di scove red (high quality) Kepler planet 
(candidate), eq.2 from lH12l is reproduced here: 



^ — ^ rt... ^ 



(7) 



j=i 



/ is evaluated for each studied stellar bin. Each planet's con- 
tribution is augmented by its geometric transit probability 
(I/Pt) to include planetary systems with non-transiting ori- 
entations as well (i.e. a detected planet with low geometric 
transit probability give a high contribution). Each planet's 
fractional contribution is calculated over the number of liigh 
photometric quality stars only (m.,j). Only those stars are 
selected from the Kepler Input Catalog (KIC) that belong 
to the analyzed stellar bin and have a high enough quality 
lightcurve where the actual planetary transit can be cer- 
tainly detected (Table [B]). A SNR value is derived from the 
edge-on planetary transit depth signa l and from the mea- 
sured scatter (ctcdpp; see eq.l in lH12l ) of the lightcurve. A 
lightcurve counts in the total number of stars if SNR>10 is 
fulfilled. 

For transiting planets, a SNR>10 transit detection (cal- 
culated from the actual transit signal depth) and orbital 
period T < 50 days are required. These planets are not 
all confirmed yet but part of the released Kepler planetary 
candidates passing an automated vetting procedure of the 
Kepler data p rocessing pipeline. They are counted as planets 
both in lH12l and in the present paper as they are assumed 



A sensitivity analysis of the WTS 17 



Stellar effective temperature bins 


3600-4100, ..., 




6600-7100K 


Stellar surface gravity, log g 


4.0-4.9 


Kepler magnitude, Kp 


< 15 


Lightcurve quality, SNR 


> 10 


Planetary radius bins, Rq 


2-4, 4-8, 8-32 


Orbital period, T 


< 50 days 


Detection threshold, SNR 


> 10 



Tabl e 6. P roperties of stellar and planetary samples considered 
in the IH12| Kepler study. 



to be actual planets with high probability (jBorucki et al.l 

To obtain stellar par ameters, the Keple r proje ct uses 
model atmospheres from ICastelli fc Kuruca (|2004l ). They 
perform a Bayesian model fitting on seven colours of 
the KIC objects determining effective temperature (Toff), 
surface gravity (logg) and metallicity (logZ) among 
other non-indepen dent model parameters simultaneously 
(JBrown et al.l 1201 Ih . Restrictive priors also ensure that pa- 
rameters remain within realistic value ranges. Temperatures 
are considered to be most reliable for Sun-like stars with dif- 
ferences from other models below 50K and u p to 200K for 
stars further away from the Sun on the CMD. lBrown et al.l 
()2011f ) call temperatures below 3750K untrustworthy. There 
are 1086 M dwarfs identified in Q2 (Table [7}. 

Using the Q2 Kepler data release, we calculate Kepler 
planetary fractions for the 0.8-10 days period range to com- 
pare with our present WTS study. We follow the steps of 
IH12| by selecting KIC objects and Kepler planets (candi- 
dates) from the Q2 data release using most criteria listed in 
Table E] but selecting planets with T < 10 days (Fig. [12]). 
We cannot filter for high quality lightc urves however as noise 
data is not readily available from lH12l . It is beyond the scope 
of this paper to process each released Kepler lightcurve and 
check its noise properties. Rather, we determine a correction 
factor for each stellar type bin. We reproduce the calculation 
for the planetary fraction in the T < 50 day case omitting 
the lightcurve SNR>10 quality criter ion and compare it to 
the published values in fig. 8 in lH12l . This reveals a correc- 
tion factor that is applied to the T < 10 day case in each 
stellar bin (Table [T]). 

For stellar bins with nonzero planets 95 per cent confi- 
dence intervals (red error bars in Fig. 1131) ar e calculated fol- 
lowing the logic of 'effective stars' from lH12l : using a Poisson 
distribution for having the actual number of planet detec- 
tions (Npi) from Npi/f 'effective stars'. Most stellar bins 
have only a few detections, so errors are heavily dominated 
by small number statistics. For bins with zero candidates (95 
per cent confidence) statistical upper limits are calculated 
in the same way as for the WTS earlier in this study. In 
this case the number of 'effective stars' are calculated as the 
integrated overall detection probability in the 0.8-10 days 
period range|j 

In Fig. [13] and Table [7] Kepler planetary fractions (/) 
are shown for short period (0.8-10 days) hot Jupiters (8- 
327?^). E.g. in the coolest stellar bin, there are 9 planets 
(2-32_R^) with period T < 50 days, counting as 320.5 occur- 



° For a Jupiter size planet, using the Kepler prior this factor is 
Pdot = 0.069. 
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Figure 12. Cumulative number of discovered Jupiter size (8- 
32R^) planets in Kepler Q2 as a function of orbital period (T) 
in different stellar temperature bins. There are no planets in the 
two coolest stellar bins (3600-4100K, 4100K-4600K). Planet oc- 
currence fractions in Fig. 1131 are calculated for the 0.8-10 days 
region, augmenting the contribution of each discovery by its geo- 
metric transit probability. 
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3.0 
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Table 7. Total number of stars (A^stars) in Kepler Q2 tempera- 
ture bins, their corresponding correction (corr.) factors (see text), 
number of Jupiter size short period planets (A^pi), their aug- 
mented contribution (A'^aug) and the occurrence ratio (/) or up- 
per limit (/95%). 



rences around 1086 M dwarfs. This y ields a fraction of 0.295 
which is the same rate published in IH12| . This means that 
the correction factor for the M dwarf bin (3600-4100K) is 1 
i.e. all 1086 M dwarf lightcurves are good enough to detect 
planets down to super-Earth sizes. As there is no Jupiter 
size planet (8-32_R0) in the 0.8-10 day period interval, we 
calculate an upper limit of the planet occurrence rate as 
3/(1086 • 0.069) = 0.04. 

We recall that in this study 2844 early (MO-2, 3400K- 
3800K) and 1679 later type (M2-4, 2960K-3400K) M dwarfs 
were identified in the 19hr field of the WTS while there 
are 1086 early M dwarfs (3600K-4100K) in the Kepler Q2 
dataset. Due to the factor of 2 higher number of early M 
dwarfs, the statistical upper limit set up by present study 
for the WTS MO-2 bin is a stricter constraint both in the 
threshold (2.9 per cent) and periodmatch approaches (3.4 
per cent) than the one that can be derived for the (early) 
Kepler M dwarfs (4 per cent, Fig. [13]). Though the WTS 
is more sensitive for Jupiters around later type M dwarfs, 
the upper limit for the M2-4 bin is slightly worse due to the 
smaller sample size; /g5% =4.3-5.2 per cent (depending on 
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Figure 13. Short period (0.8-10 day) liot Jupiter planetary oc- 
currence fractions (x) and upper limits (4,) in case of null detec- 
tions in the WTS (gree n), Kepler Q2 ||H12|) (red) and from the 
HARPS planet search l lBonfils et ai]|201ll) (blue). WTS upper 
limits are shown for the MO-2, M2-4 bins and also for the whole 
MO-4 sample (green arrows). Vertical error bars on WTS upper 
limit markers cover the uncertainty from the threshold and pe- 
riodmatch interpretations (Table [5}. Horizontal bars (I — \) show 
temperature bin widths, horizontal arrows (f-^) mark estimated 
uncertainty of bin edges where available. All upper limits and 
error bars of nonzero fractions are for 95 per cent confidence. 



the period prior), while for the overall MO-4 sample, it is 
1.7-2.0 per cent. 

The WTS upper limits sit above the measured Hot 
Jupiter fractions found for hotter stars by Kepler which have 
much less uncertainty from much larger samples (bins with 
higher than 4600K in Fig. [13]). With our analysis of the 
WTS, we cannot rule out a similar Hot Jupiter occurrence 
rate around M dwarfs (MO-4). Our main conclusion is that 
there is currently no evidence to support the argument that 
Hot Jupiters are less common around M dwarfs than around 
solar-type stars. 

To push our constraints to lower levels we will need to 
complete the remaining three fields of the WTS, or wait 
for the completion of oth er survey s such as Pan-P lanets 
(JKoppenhoefer et al.llJOOgl ) or PTF (JLaw et al.ll2012l ). 

The completed WTS survey contains a total of almost 
15,000 MO-4 dwarfs. A confirmed null detection for hot 
Jupiters in the larger sample would push our upper limit 
on the HJ fraction down by a factor 3.2, i.e. to /95% =0.5- 
0.6 per cent for MO-4 stars. This will then overlap with the 
values found for hotter stars in Fig. 1131 

We note tha t the recent Kepler d iscovery and confirma- 
tion of KOI-254 ([Johnson et al.ll2012l ). a hot Jupiter planet 
around an M dwa rf is beyond the parameter space con- 
sidered in the IH12| study thus cannot be directly incorpo- 
rated into our comparison ( objec t is fainter). Nevertheless, 
if this planet were part of the lH12l sample, this sole detection 
would mean a l.llii per cent occurrence rate which is lower 
than the upper limits b ut still compa t ible w ith the findings 
of the present and the iBonfils et al.l (|201ll ) RV study. We 
note that this is an overestimation of the real weight of this 



planet detection as there are more cool hosts in KIC down 
to Kp =15.979 but linked to unknown, possibly lower survey 
sensitivity as well. 

There are some caveats in our comparison, (i) The 
early M dwarf (MO-2) temperature bins are not exactly 
th e sam e in the present WTS study (3400K-3800K) and 
in |H12| (3600K-4100K) thus the WTS bins are apparently 
~ 250K cooler. However, we cannot use the same method 
to re-estimate temperatures for Kepler sources. Hopefully, 
as more follow-up data is released this will become possible. 
(ii) The sample of Kepler objects used to calculate plane- 
tary occurrences is restricted to bright objects [Kp < 15) 
which is different from the magnitude range of WTS objects 
studied here (J < 17). It is possible that the Kepler early M 
dwarf sample is from a different popula tion than the sample 
in the present study (jMann et al.ll2012l ). (iii) The sensitivity 
for Jupiters in the WTS was determined for a conservative 
10-Re radius while hot Jupiters may have larger radii. The 
Kepler study uses the 8-32_R® radius range for hot Jupiters. 



7.1 Comparison with RV studies 

iBonfils et al.l (|201ll ) analyze RV data observed by the 
ESO/HARPS spectrograph for planet signatures in their M 
dwarf (MO-6) sample of 102 nearby (<llpc, V < 14) stars 
and find no Hot Jupiters. According to their sensitivity con- 
siderations, this sample is equivalent to 96.83 effective stars. 
Using our 95 per cent confidence level, their null detection 
of hot Jupiters (1-10 days) implies a 3.1 per cent upper limit 
on the planetary occurrence rate around M dwarfs, in good 
agr eement with o u r WT S results. 

IWright et all (|2012l ) find that 1.2±0.38 per cent of 
nearby Solar-type stars host hot Jupiters. The ir results are 
in good agreement with p r evious RV works (JMarcv et al.l 
l2005l : ICumming et al.lboOSi : iMavor et al.ll20lil ). Occurrence 
rates found by tr a nsit studies are systematically lower 
(|Gould et all (|2006| ). IH12| ). around 0.5 per cent. 

As the rates around G dwarfs are about the same (RV) 
or lower (transit) than the upper limits reached by the WTS 
around M dwarfs in present study, we cannot rule out planet 
formation scenarios in this paper. We admit that this is a 
very brief comparison only, avoiding the discussion of the 
methods and error levels used by the studies cited above. 



8 SUMMARY 

The WFCAM Transit Survey has observed 950 epochs for 
one of its four target fields. These data were collected over 
more than three years and will provide a valuable resource 
for general studies of the photometric and astrometric prop- 
erties of large numbers of objects in the near-infrared. In this 
paper we determined the sensitivity of this dataset to short 
period (<10 day) Jupiter and Neptune sized planets around 
host stars of spectral type M0-M4. We identify and classify 
two subsamples of M dwarfs: M0-M2 comprising 2844 stars, 
and M2-M4 comprising 1679 stars. We compare Dartmouth 
and NextGen models to derive estimates of Tcb for the sam- 
ple, and demonstrate that reddening effects are small for 
the 19-hour field. This forms one of the largest samples of 
M dwarfs targeted by any dedicated transit survey, and cur- 
rently the only one working in the near-infrared. 
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We have described how our multi epoch WFCAM data 
are filtered and processed to produce lightcurves with me- 
dian RA4S ~ 4 mmag down to J = 14 and RMS ~ 17.5 
mmag at J = 17. We found tlie RMS for the briglitest stars 
suffers a significant contribution from systematic noise at 
the Js 3 mmag level. The origin of this is not determined, 
but we note that we see similar (albeit smaller) effects in the 
optical (e.g. Ilrwin et al.ll2007l ). We leave for future work an 
investigation of the flatfield and near-infrared background 
as two potential contributors to the systematics. 

We performed Monte-Carlo simulations on the WTS 
lightcurves, injecting and recovering fake transit events. We 
generalise each of the stellar and planet populations into 
two distinct radius regimes, giving us 4 scenarios for tran- 
sit depths: Jupiters and Neptunes, around M0-M2 stars and 
M2-M4 stars. We investigate the resultant signal-to-noise of 
the recovered events and compare this to the survey detec- 
tion thresholds. We also investigate our sensitivity to these 
4 scenarios as a function of orbital period and stellar mag- 
nitude. Our analysis of the simulations enables us to place 
constraints on the incidence of hot Jupiters. 

With 95 per cent confidence and for periods <10 days, 
we showed that fewer than 4.3-5.2 per cent of M2-M4 dwarfs 
host hot Jupiters. Constraints are even stronger for earlier 
spectral types, and fewer than 2.9-3.4 per cent of M0-M2 
dwarfs host hot Jupiters, while for the overall MO-4 sample, 
it is 1.7-2.0 per cent. An analysis of the Kepler Q2 data 
shows that the WTS provides more rigorous upper limits 
around cooler objects, thanks to the larger size of our sample 
(2844 M dwarfs in the WTS compared to 1086 in the Kepler 
sample under consideration). We compare these upper limits 
to the measured hot Jupiter fraction around more massive 
host stars, and find that they are consistent, i.e. we cannot 
rule out similar planet formation scenarios around at least 
the earlier M dwarfs (MO-4) at the moment. 

We used our simulations to demonstrate that the WTS 
lightcurves are sensitive to transits induced by hot Neptunes 
in favorable cases, but that the ability to distinguish them 
from false alarms (astrophysical or systematic) is limited by 
a number of factors. Firstly the transit events have signifi- 
cantly lower to signal-to-noise than events arising from hot 
Jupiters, and thus the contamination by false alarms (from 
systematic noise) is higher at the required detection thresh- 
olds. Secondly, the recovered periods for the simulated hot 
Neptunes are dominated by spikes at multiples of 1 day, not 
seen so strongly for the larger planets. Without reliable peri- 
ods, any attempts at follow-up photometry and spectroscopy 
for detected candidates will be hardly viable, especially for 
such faint target stars. 

The data presented in this paper represent one quarter 
of the planned WFCAM Transit Survey (the other 3 fields 
currently lack sufficient coverage to include in the analysis). 
If the WTS is completed then we expect our statistical con- 
straints to improve by a factor of 3-4 (dependent on the 
actual numbers of M dwarfs in the remaining three fields). 
If no hot Jupiters are found, our upper limits get to the 0.5 
per cent level for M0-M4 dwarfs, and may ultimately reach 
a significantly lower level for hot Jupiter occurrence than is 
measured around G dwarfs by other surveys. 
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